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ABSTRACT 

We have obtained near-infrared imaging of 58 galaxy groups, in the redshift range 0.1 < 
z < 0.6, from the William Herschel Telescope and from the Spitzer IRAC data archive. 
The groups are selected from the CNOC2 redshift survey, with additional spectroscopy from 
the Baade telescope (Magellan). Our group samples are statistically complete to -KVega = 
17.7 (INGRID) and [3.6^to]ab = 19.9 (IRAC). From these data we construct near-infrared 
luminosity functions, for groups in bins of velocity dispersion, up to 800 km s _1 , and redshift. 
The total amount of near-infrared luminosity per group is compared with the dynamical mass, 
estimated from the velocity dispersion, to compute the mass-to-light ratio, M200 /Lk- We find 
that the M200 /Lk values in these groups are in good agreement with those of their statistical 
descendants at z = 0, with no evidence for evolution beyond that expected for a passively 
evolving population. There is a trend of M200/ Lk with group mass, which increases from 
M 2 oo/Lk~ 10 for groups with a < 250 km s _1 to M 2 oo/Lk~ 100 for 425 km s _1 < a < 
800 km s . This trend is weaker, but still present, if we estimate the total mass from weak 
lensing measurements. In terms of stellar mass, stars make up > 2 per cent of the mass in the 
smallest groups, and < 1 per cent in the most massive groups. We also use the near-infrared 
data to consider the correlations between stellar populations and stellar masses, for group and 
field galaxies at 0.1 < z < 0.6. We find that fewer group galaxies show strong [Oil] emission 
compared with field galaxies of the same stellar mass and at the same redshift. We conclude 
that most of the stellar mass in these groups was already in place by z ~ 0.4, with little 
environment-driven evolution to the present day. 
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1 INTRODUCTION 



In hierarchical models of structure formation, small groups of 
galaxies represent an important environment that connects active, 
star-forming field galaxies to quiescent cluster galaxies. In this 
context they may be the environment in which the star formation 
histories of galaxies are dramatically altered. Since a large frac- 
tion of galaxies in the nearby Universe exist within groups (e.g. 
iHuchra & GelleJl982tlEke et all20043) . such group-driven trans- 
formations could have a dominant influence on galaxy evolution, at 
least since z ~ 1. However, galaxy groups are not well understood 
relative to field galaxies and larger clusters, especially beyond the 
local universe. Their low-contrast relative to the background means 
that deep, highly complete redshift surveys are needed to compile 



large enough samples of member galaxies, and to minimize the ef- 
fects of selection biases. 

Now that the global history of star formation over the past few 
billion years is well estab lished (e.g. lHeavens et all2004llHopkinsl 
12004 Iluneauetalj l2005). the next important goal is to trace star 
formation and efficiency in different clustering environments (e.g. 
iBaloeh et alll2004 lEke e7aTll2004rJ : Icooper et al Jl200rj) . An im- 
portant integral quantity that traces this evolution is the relative 
amount of luminous matter in groups of different masses. Stud- 
ies at optical wavelengths have shown that, locally at least, the 
ratio between dynamical mass and optical light (M2oo/L opt ) in- 
creases strongly with mass, up to systems w ith masses of a few 
times 10 13 M(T), typical of poor clusters (e.g. iGirardi et alj|2002t 
Eke et al. 2004b). Semi-analytic models of galaxy formation with 
reasonable feedback schemes are in good agreement with this trend 
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lEke et aljE004bl) . Recently. IFarker et alJ feOOSt) have made mea- 
surements of M200 /iopt in galaxy groups selected from the second 
Canadian Network for Observational Cosmology (CNOC2) field 
galaxy redshift survey ( Carlberg et al JI2001I). with masses deter- 
mined by weak lensing (Hoek stra et alJl200ll) . They, too, find a 
trend of increasing Mzoo/Lopt that is consistent with that mea- 
sured at lower redshift. 

A difficulty in interpreting the above results is the sensitivity 
of optical luminosity to recent star formation , which has a strong 
environmental dependen ce (e.g. iBalogh et alj|2004l iBlanton et all 
l2005tlBaldrv et all20()3) . Near-infrared (NIR) luminosities are su- 
perior in this respect, because they trace stellar mass, indepen- 
dent of star formatio n history, to within a factor of a few (e.g. 
Bell & de Jonj200lh . A few near-infrared studies of nearby groups 
have been undertaken fcin et aljEbol iRamella et all feoCM). and 
these observations generally find an increase in M200 / Lk (the ra- 
tio between dynamical-mass and near-infrared light) with mass, 
comparable to that seen at optical wavelengths. However, until now 
such measurements have not yet been reported for more distant 
groups. 

It is the purpose of this paper to measure M200 / Lk for groups 
at 0.1 < z < 0.6 for the first time, using both ground-based 
(K-band) and Spitzer (3.6^im) imaging of 58 galaxy groups from 
the CNOC2 survey. The paper is outlined as follows: In Sj2|below 
we review the CNOC2 survey, discuss our observations, data re- 
duction techniques and methods of group member selection. We 
present the luminosity functions and M200/LK as a function of 
group redshift and velocity dispersion in ij3| The implications of 
these measurements, interpreted as stellar masses, are discussed in 
Sj4] and our conclusions are summarized in §5\ Throughout this 
paper we assume a cosmology with matter density Q m = 0.3, 
energy density = 0.7, and current Hubble constant H a = 
lOO/i km s- 1 Mpc" 1 with h = 0.75. 



2 THE DATA 

2.1 The redshift surveys 

Our galaxy group sample is derived from the CNOC2 redshift sur- 
vey, which was conducted at the Canada-France-Hawaii telescope 
using the multiobject spectrograph (MOS), with the pri mary goal 
of studying galaxy populations, clustering, and evolution lYee et all 
2000). This spectroscopic and photometric survey consists of ap- 
proximately 6000 galaxies with a measured redshift, distributed in 
four widely separated patches covering a total of ~ 1.5 deg 2 of the 
sky. Although the survey area is small relative to more recent red- 
shift surveys, the use of four noncontiguous patches ensures that 
our group sample is not likely to be strongly biased due to cosmic 
variance. The survey includes U, B, V, R c , I photometry, and the 
initial spectroscopic sample was selected from the R c (hereafter 
just R) images. The colours are k-corre cted to the rest-frame using 
the template-fitting code of lBlanton & Roweisl j2006). Each MOS 
pointing was observed with two masks, to increase completeness in 
dense fields. There are redshifts for about 50 per cent of galaxies to 
R = 21.5; this completeness is magnitude-dependent, as slits were 
primarily allocated to brighter galaxies. 

Recently we re-observed 20 fields in three of the four pa tche s 
with the LDSS2 spectrograph on Magellan (Wilman et al. 2005a). 
There were two main aims of this additional spectroscopy. The first 
was to cover each group with sufficient masks (between one and 
three) to achieve near- 100 per cent completeness at bright magni- 
tudes. This is important because measurements of quantities like 



total stellar mass will be strongly affected by the brightest galaxies, 
and their small number per group will make statistical corrections 
unreliable on an individual group basis. Furthermore, the use of a 
larger telescope and better spectrograph allows us to obtain red- 
shifts for fainter galaxies, and we obtain a good statistical sampling 
of galaxies as faint as R = 22. The main advantage for our present 
purposes is that this increases the number of members per group, 
allowing us to obtain better measurements of the group geometry 
and velocity dispersion. 

From the CNOC2 and Magellan spectra we also measure the 
rest-frame equivalent width of the [Oil] emission line, ^([Oll]), 
as an indicator of ongoing star formati on. This line str ength is mea- 
sured following the definition of Bal ogh et alJ ll997l) . In the fields 
covered by our LDSS2 spe ctroscopy, the line me asurements are the 
same as those presented bv lWilman et aljl200 5a^: for the rest of the 
CNOC2 sample, the measurements will be presented in a future pa- 
per (Whitaker & Morris, in preparation). 

2.2 Galaxy Group Membership 

Virialized galaxy groups within the CNOC2 patches were initially 
identified in redshift space by Carlberg et al. 1 200 1 ) using an iter- 
ative friends-of-friends algorithm on the catalogued galaxies. Over 
200 groups were found, with an average o f 3.8 co nfirmed members 
per group. The search by Carl berg et alJ 1200 D limited the sam- 
ple's fc-corrected and evolution-compensated absolute luminosity 
to brighter than Afjj e = —18.5 mag. To take advantage of the in- 
creased depth and completeness afforded us by the Magellan data, 
in this paper we redefine the original group membership, in a man- 
ner similar to that described in Wilman et al. (2005a). To summa- 
rize briefly, an iterative procedure is implemented to select galaxies 
with redshifts within twice the group velocity dispersion and with 
a transverse distance from the group centre which is within 1/5 
of the line-of-sight distance. In each iteration, the velocity disper- 
sion is recomputed using the Gapper estimate iBeers et al]ll990h . 
and the centre is recomputed as the luminosity-weighted geomet- 
ric centre of the group members. This procedure is applied to the 
whole group sample, including those that were not re-observed with 
LDSS2, so that membership is consistently defined. 

To compute the virial mass of all the galaxy groups we used a 
theoretical estimate based on the velocity dispersion, a. Assuming 
that the groups are in dynamical equilibrium, the total mass con- 
tained within a volume of |7rJ?200 is given by 



3 2 

7M200 = -^RlOQO , 
(jr 



where the "virial radius" -R200 is 



Hi 



V3a 



(1) 



(2) 



10H o (1 + z) 1 ' 5 ' 
Although we are well aware that velocity dispersion is not a per- 
fect indicator of mass, especially for these poor systems with few 
members, we keep this definition so that w e can di rectly and easily 
compare our results with the literature (e.g. Ramella et al. 2004). In 
Section l4~T1 we will consider how the use of available weak-lensing 
masses affects our results. 

As described in the following sections, we have obtained near- 
infrared imaging for 58 groups, and these are the only groups we 
will consider in the remainder of this paper. The list of groups and 
their properties are given in Table Q for those groups with only 
data from the original CNOC2 survey, and in Table|2|for those re- 
observed with Magellan. Included in the tables are the group coor- 
dinates, redshift, intrinsic velocity dispersion (i.e., accounting for 
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Table 1. Properties of galaxy grou ps from the origina l ICarlberg et all jJlOOl) sample (i.e. not reobserved at Magellan) with good near-infrared coverage. 
Column 1 gives the group id from Carlberg et al. 1 2001), modulus factors of 100 added to distinguish groups in different patches. Columns 2-4 give the central 
position and redshift of the group. Column 5 is the estimated intrinsic velocity dispersion, and column 6 is the number of galaxies used to obtain this number. 
Columns 7 and 8 are the virial radius and mass calculated from the velocity dispersion. Entries marked with an * are upper limits. The numbers in columns 
(9), (10) and (11) are the number of galaxies with redshifts, within i?200 in the R— band, Spitzer and INGRID catalogues, respectively. 



(1) 


(2) 


(3) 


(4) 


(5) (6) 


(7) 


(8) 


(9) 


(10) 


(ID 


Group 


a 


S 


z 


c AT CT 


-R200 


M200 


N z 


N z 






(J2000) 


(J2000) 




(km s -1 ) 


(kpc) 


(10 12 M e ) 


(R) 


(Spitzer) 


(INGRID) 



1 


222.4258728 


9.050009727 


0.1647 


211 ± 95 


8 


388 ± 174 


1.21 ± 1.63 


5 


5 


4 


9 


222.2241516 


8.945289612 


0.2616 


194 ± 61 


8 


316 ± 99 


0.83 ± 0.78 


4 


3 


3 


1 1 


222.2881775 


8.830399513 


0.2708 


185 ± 46 


12 


298 ± 75 


0.71 ± 0.54 


4 


4 





13 


222.3945923 


8.905659676 


0.2712 


397 ± 68 


1 1 


640 ± 110 


7.05 ± 3.62 


4 


4 





15 


222.1999817 


8.963080406 


0.3071 


308 ± 219 


4 


476 ± 338 


3.14 ± 6.71 


3 


2 


3 


16 


222.5902405 


9.103640556 


0.3063 


254 ± 64 


5 


392 ± 100 


1.76 ± 1.34 


4 


4 





17 


222.097641 


8.784229279 


0.3086 


395 ± 92 


8 


610 ± 141 


6.66 ± 4.62 


3 





3 


18 


222.0555878 


8.9292202 


0.3234 


261 ± 264 


4 


397 ± 401 


1.89 ± 5.74 


4 





4 


19 


222.5530853 


8.960080147 


0.3248 


330 ±115 


6 


499 ± 174 


3.79 ± 3.96 


4 


4 





21 


222.4801636 


9.646129608 


0.3483 


204 ± 73 


7 


301 ± 107 


0.87 ± 0.93 


6 





5 


29 


222.4464874 


8.852370262 


0.3736 


309 ± 158 


7 


443 ± 227 


2.95 ± 4.53 


4 


4 





30 


222.4998779 


8.82020092 


0.3938 


335 ± 87 


8 


469 ± 122 


3.67 ± 2.85 


5 


4 





36 


222.3757324 


9.153150558 


0.4701 


265 ± 308 


4 


343 ± 399 


1.68 ± 5.85 


4 


4 





202 


36.58484268 


0.4656900167 


0.1883 


341 ± 288 


5 


607 ± 513 


4.91 ± 12.45 


4 


4 





206 


36.37388992 


0.2616599798 


0.2284 


261 ± 63 


1 1 


442 ± 106 


2.09 ± 1.51 


6 


4 





208 


36.48381424 


0.1982000023 


0.2686 


592 ± 165 


8 


956 ± 267 


23.38 ± 19.54 


7 


7 





217 


36.44522095 


0.3269200325 


0.3082 


696 ± 276 


5 


1074 ± 425 


36.31 ± 43.12 


5 


5 





221 


36.52441406 


0.07746999711 


0.3578 


73 ± 95 


6 


107 ± 139 


0.04 ± 0.16 


1 


1 





229 


36.69626617 


0.142960012 


0.3833 


540 ± 139 


13 


766 ± 197 


15.56 ± 12.02 


3 


3 





234 


36.5413208 


0.4848300517 


0.3974 


284 ± 111 


7 


397 ± 155 


2.24 ± 2.62 


3 


3 





239 


36.36482239 


0.2544400096 


0.4083 


608 ± 96 


6 


841 ± 133 


21.72 ± 10.29 


6 


6 





301 


141.086853 


36.9005127 


0.1072 


66 ± 86 


6 


130 ± 170 


0.04 ± 0.15 


2 


2 





302 


140.2914429 


36.71305847 


0.112 


123* 


3 


242* 


0.26* 


3 





3 


303 


141.0169373 


37.0852623 


0.1913 


331 ± 98 


7 


587 ± 175 


4.48 ± 4 


8 


6 





304 


140.8880463 


37.38243866 


0.1914 


121 ± 159 


3 


214 ± 282 


0.22 ± 0.86 


3 





2 


305 


140.9043732 


36.91950989 


0.2025 


243 ± 104 


5 


425 ± 182 


1.75 ± 2.24 


3 


3 





313 


140.9319 


37.01506042 


0.2335 


232 ± 97 


7 


391 ± 164 


1.47 ± 1.85 


4 


4 





315 


141.0404053 


36.99619293 


0.2435 


243 ± 45 


13 


404 ± 74 


1.66 ±0.92 


5 


5 





317 


140.9354401 


37.34311676 


0.2449 


280 ± 66 


12 


465 ± 109 


2.54 ± 1.79 


9 





6 


329 


140.9833832 


37.25820923 


0.3223 


487 ± 79 


4 


739 ± 120 


12.22 ±5.94 


1 





1 


333 


140.8763885 


37.14134979 


0.3218 


455 ± 143 


9 


692 ±217 


10.02 ±9.42 


6 


5 


4 


336 


140.8822021 


36.8615303 


0.3629 


463 ± 531 


4 


672 ± 771 


10.05 ±34.58 


4 


4 


4 


337 


140.9049683 


37.7922287 


0.3729 


624 ± 184 


7 


895 ± 264 


24.31 ±21.54 


6 





4 


344 


141.0652008 


36.88687134 


0.3735 


277 ± 77 


12 


397 ± 110 


2.12 ± 1.77 


3 


3 





346 


141.004715 


36.78276443 


0.3733 


393 ± 30 


26 


563 ± 44 


6.06 ± 1.41 


12 


12 





348 


140.3150482 


36.71141052 


0.3789 


82 ± 106 


4 


117 ± 151 


0.05 ±0.21 


1 





1 


349 


140.6838989 


36.95294952 


0.38 


498 ± 108 


8 


710 ± 154 


12.3 ±8.02 


3 





3 


350 


140.5788269 


36.81409073 


0.3783 


285 ± 85 


7 


407 ± 122 


2.31 ±2.08 


3 





3 


351 


140.965744 


37.26145172 


0.3795 


143 ± 88 


6 


203 ± 126 


0.29 ±0.53 


2 





2 


352 


141.0719757 


36.94703293 


0.3791 


386 ± 166 


4 


550 ± 236 


5.72 ± 7.37 


3 


3 





355 


140.6422424 


36.72232056 


0.3907 


349 ± 220 


5 


491 ±310 


4.18 ± 7.9 


3 





2 


357 


140.3446808 


36.62179947 


0.3909 


133 ± 107 


7 


187 ± 151 


0.23 ±0.56 


2 





2 


358 


140.7977142 


37.12773132 


0.3908 


287 ± 53 


11 


404 ± 75 


2.32 ± 1.29 


6 





6 


359 


140.9085388 


37.121418 


0.3911 


266 ± 127 


4 


375 ± 179 


1.85 ± 2.65 


3 


3 


3 


361 


140.7460632 


36.83348846 


0.4275 


118 ± 156 


3 


159 ±211 


0.15 ±0.61 


2 





2 


365 


141.1464691 


37.03060913 


0.473 


122 ± 163 


3 


158 ±211 


0.16 ±0.66 


3 


3 





366 


140.6099243 


36.71389771 


0.4728 


455 ± 191 


9 


588 ± 246 


8.49 ± 10.68 


4 





4 



velocity uncertainties), virial mass, and the number of group mem- 
bers we identify within -R200 above the R— completeness limit (see 
i|2.3t . The velocity dispersions were computed from all galaxies 
within 1/h Mpc of the group centre (unless there were less than 
four members, in which case all galaxies were used); the number of 
galaxies contributing to the measurement of a is given in the table 
as TVtj, as is the the number of group members with either INGRID 
or Spitzer observations, as described in Sections l2.4l and l2~5l 



2.3 Redshift completeness 

The CNOC2 redshift survey is ~ 45 per cent complete on aver- 
age, with a depe ndence on R magnitude. For the original survey, 
lYee et all J2000) computed statistical weights that correct for this 
incompleteness, and are valid to a limit of R = 21.5. Th e addi- 
tional, Magellan spectroscopy we obtained <Wilmanetail E005a) 
improves the depth around selected galaxy groups to R = 22. It 
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Table 2. As TableQ but for the groups reobserved at Magellan <Wilman et all2005 J) , with good NIR data. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Group 


a 


5 


z 


a 




A200 


M200 


N z 


N z 


N z 




(J2000) 


(J2000) 




(km s -1 ) 




(kpc) 


(1O 12 A4 ) 


(R) 


(Spitzer) 


(INGRID) 


23 


222.3715668 


9.511750221 


0.3515 


510 ± 113 


8 


749 ± 167 


13.58 ±9.06 


6 





6 


24 


222.2642517 


9.116889954 


0.3592 


121* 


10 


177* 


0.18* 


6 


5 





27 


222.4210205 


9.037179947 


0.3725 


181 ± 148 


3 


260 ±212 


0.59 ± 1.45 


3 


2 


3 


28 


222.5953522 


9.018850327 


0.3728 


161 ± 79 


6 


231 ± 113 


0.42 ±0.61 


4 


4 


4 


31 


222.3106995 


9.188610077 


0.3929 


450 ± 455 


4 


632 ± 640 


8.94 ± 27.14 


3 


2 





32 


222.4871063 


8.929120064 


0.3939 


591 ± 120 


6 


829 ± 169 


20.22 ± 12.32 


5 


5 





33 


222.3906708 


9.480049133 


0.4066 


125 ±81 


6 


173 ± 113 


0.19 ±0.37 


2 





2 


37 


222.3853607 


9.073459625 


0.4713 


236 ± 79 


11 


305 ± 102 


1.18 ± 1.19 


4 


4 


4 


38 


222.3486328 


8.980949402 


0.5106 


773 ± 77 


14 


962 ± 95 


40.17 ±11.93 


10 


8 





39 


222.3684692 


9.496768951 


0.5366 


461 ± 104 


13 


559 ± 126 


8.29 ±5.59 


8 





8 


227 


36.62685013 


0.2123200148 


0.3635 


341 ± 144 


7 


495 ± 209 


4.02 ± 5.08 


4 


4 






is difficult to combine the two samples in a statistically fair way; 
we thus keep them separate, and will refer to them as the R21.5 
and R22 samples, respectively. These samples are exclusive: that 
is, galaxies and groups in the R22 sample are removed from the 
R21.5 sample, even though, of course, they formed part of the origi- 
nal CNOC2 survey. The statistical weights previously computed for 
both these surveys are magnitude- and position-dependent, but are 
computed from ensemble averages that may not appropriately ac- 
count for any variation in completeness that depends on group rich- 
ness, since groups with more members may be more incompletely 
sampled, despite the multiple-mask observing strategy. We have 
therefore recomputed the magnitude weights for group members as 
follows. First, the group sample is divided by velocity dispersion 
into those with a < 200 km s -1 , 200km s _1 < a < 500km s -1 
and a > 500 km s~ 1 . For each of these three subsamples, we com- 
pute the fraction of galaxies within -R200 that have a redshift, as a 
fu nction of R m agnitude only, following the procedure described 
in lWilman et alJ l2005al) . but omitting the radial weight. For field 
galaxies, we use a simple magnitude-dependent weight that is com- 
puted independently for the R21.5 and R22 samples. Despite all 
these efforts, our results are only weakly dependent on the details 
of the weighting scheme. In particular, we find retrospectively that 
the incompleteness in the most massive groups is at most ~ 10 
per cent greater than the average. Magnitude weights are generally 
less than 1.5 for galaxies brighter than R=20, and increase up to a 
maximum of 5 (i.e. 20 per cent completeness) at R=22. 



2.4 K-band Observations and Data Reduction 

Over the course of three nights during March 9-13, 2001, 25 fields 
centred on groups in the 9h and 14h CNOC2 patches were observed 
with the Isaac Newton Group Red Imaging Device (INGRID) on 
the 4.2m William Herschel Telescope (WHT). Each field was ob- 
served using a nine-point dither pattern with exposures of 10— 15s 
in K 3 . The total integration for each field varied, but was typ- 
ically about 12 minutes. The data were re duced using the IPIPE 
NIR reduction pipeline iGilbank et ail20 03). Briefly, this involved 
bad pixel-masking, flat fielding and sky-subtraction using a local 
flat field, and offset-finding and stacking of the sky-subtracted im- 
ages to create a first- pass mosaic. An object mask was generated 
using SEXTRACTOR iBertin & Arnoutslll99fj) to detect objects in 
this mosaic, and the sky-subtraction process was repeated after ap- 
plying this mask to the individual exposures. This two step sky- 
subtraction technique avoids over-subtraction of the sky, as faint 



objects undetected in a single frame can otherwise lead to an over- 
estimate of the sky level. The second pass sky-subtracted frames 
were then combined into a final mosaic using a 3-a clipped mean. 

Our magnitude zeropoints (on the Vega system) and astromet- 
ric solution were calibrated using the Two Micron All Sky Survey 
(2MASS; see ljarrett et aljEo OO) point-source catalog i^-band 4" 
standard aperture magnitudes. Although the extended source cata- 
log would have been a preferable source of photometry, there were 
not enough matching objects in our fields to obtain a reliable cali- 
bration. We used SEXTRACTOR v2. 3. 2 to obtain 4" aperture mag- 
nitudes from our observations, and these were then compared to the 
matching 2MASS point-source catalog object magnitudes. For ev- 
ery field, the median of the differences between the SEXTRACTOR 
and 2MASS magnitudes was taken to be the zeropoint shift 1 . 

For our analysis we adopt th e MAG_BEST photometry output 
from SEXTRACTOR. This gives a Kronl (19801) magnitude if there 
is no danger of a nearby object biasing the magnitude measurement 
by more than 10 percent. If such crowding is an issue, a corrected 
isophotal magnitude is used instead. The perimeter of all images 
were trimmed such that most of the underexposed regions were not 
included in our analysis. We discarded all detected objects with 
SEXTRACTOR flags greater than 8. Two of the 25 images obtained 
were rejected entirely; one due to poor image quality, and the other 
due to poor seeing (~ 1.7"). Obj ects were then matched with the 
CNOC2 photometric catalogues <Yeeetalj f200C). Unfortunately, 
the global astrometry of the optical catalogues is not better than 
~ 5". This effect was mitigated by adjusting the centroid of each 
field independently; over the relatively small INGRID field-of-view 
(4.1') the astrometry was good enough to allow reliable matching. 

In the top panel of FigureQ we show the K-R colour as a func- 
tion of R magnitude, for all objects classified as galaxies. From the 
absence of blue galaxies fainter than R ~ 20 we determine that 
the photometric magnitude limit is approximately K = 19.4. We 
will show below CS I2-6I that this is much fainter than the complete- 
ness of the spectroscopic sample (shown as the dotted line), so all 
galaxies with redshifts are well-detected in the near-infrared. 



1 In practice this was done using the brightest sources, but fainter objects 
were included if there were enough of them to overcome the larger uncer- 
tainties on their individual magnitudes. 
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Figure 1. The optical-NIR colours of galaxies in the INGRID sample (top 
panel) and (for clarity) one-third of the Spitzer sample (bottom panel), as 
a function of R magnitude. Crosses represent galaxies with spectroscopic 
redshifts. The solid lines represent the photometric completeness of the IR 
data, which are K = 19.4 and SI = 23.4, respectively. The dotted lines 
represent the limits of K = 17.7 and 51 = 19.9 that correspond to the 
spectroscopic completeness limit for z < 0.6 (see § 12.61 . 



2.5 Spitzer observations 

Data taken with the Spitzer telescope Infrared Array Camera 
(IRAC) were obtained from the archive (GO program 64, PI G. 
Fazio). The pipeline-reduced images were used, and objects were 
detected in IRAC band 1 (AB magnitudes at 3.6/im, hereafter de- 
noted SI) with SEXTRACTOR v2.4.4. The detection criterion was 
a minimum of four pixels above a threshold of 1.5 times the global, 
rms background, with no filtering. The deblending algorithm used 
a minimum contrast of 0.00002, with 64 sub-thresholds. These pa- 
rameters were chosen to optimize the deblending for matching with 
the optical catalogues. The object fluxes were determined from the 
improved MAG-AUTO magnitude, which is recommended for the 
most recent versions of SEXTRACTOR. 

Again, the relatively poor astrometry of the optical catalogues 
causes some difficulty in cross-correlation. For the 14h patch, we 
have improved the astrometric solution based on several Hubble 
Space Telescope pointings. However, in the 2h and 9h fields, the 
astrometric solution varies by ~ 5" across the overlapping IRAC 
field-of-view. Thus a larger matching radius was chosen (4" in the 
2h field and 5" in the 9h field), which causes an estimated ~ 10 per 
cent of sources to be incorrectly matched. The photometric com- 
pleteness is approximately S 1=23.4, as shown in the bottom panel 
of Figure Q As is the case with the K-band data, this is several 
magnitudes below our spectroscopic limit. 



2.6 Completeness 

We need to define the magnitude limits in both K and 51 such that 
the spectroscopic sample of galaxies above these limits is statisti- 
cally complete. This will be the case if R < 21.5 (for the R21.5 
sample) or 7? < 22 (for the R22 sample) for all galaxies brighter 
than these limits. Figure|2|shows the optical-NIR colours as a func- 



Figure 2. Observed optical-IR colours are shown as a function of redshift, 
for all galaxies with spectroscopic redshifts. Filled symbols represent group 
members. The R band magnitudes come from the CNOC2 catalogue, while 
the S 1 magnitudes (bottom panel) are from the Spitzer IRAC data, and the K 
magnitudes (top panel) are from INGRID. The lines are Bruzu al & Charlo< 
(2003) models for a solar metallicity, 13.7 Gyr model. The dotted line is 
dust-free, while the solid line includes t v = 1 magnitude of extinction, and 
thus provides a good estimate for the maximum possible colour for most 
normal galaxies. 



tion of redshift, for the INGRID and Spitzer samples. The solid line 
shows a model of a solar metallicity, d ust-reddened, 13.7 Gyr old 
stellar population from the lBruzual & CharlotH200 3) models. Here 
an d throughout the pape r, we use the two-component dust model 
of ICharlot & Falll 12000). assuming a visible optical depth to the 
youngest stars (with age < 10 7 yr) of r„ = 1, and 30 per cent 
of this value toward older stars. The model provides a reasonable 
match to the red envelope of the data. The data and model indicate 
that, at z < 0.6, all normal galaxies should have (R — SI) ^ 2.1 
and (R — K) < 4.3. Galaxies redder than this (of which there 
are many in our photometric sample) are almost certainly at higher 
redshift. Given the R-completeness limits of the R21.5 and R22 
samples, we therefore have the following NIR limits: 



Knm = 



Sllin 



17.2 #21.5 

17.7 R22 

19.4 #21.5 

19.9 #22 



(3) 



(4) 



In all cases, these magnitudes are well above the photometric com- 
pleteness limits. Therefore, applying the weights appropriate to the 
R21.5 or R22 sample will provide statistically complete samples to 
these limits. 

To compute the stellar luminosity of all group members, we 
need finally to account for the fact that some groups are only par- 
tially covered by Spitzer or INGRID imaging. For Spitzer, this gen- 
erally occurs for groups that lie near the edge of the IRAC image. 
The INGRID images were centred on particular groups, but could 
still be incomplete if the virial radius (Equation|2j is larger than the 
field of view, or if there are background (or foreground) groups in 
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Figure 3. Model k-corrections, from a suite of Bruzual & Charlo] I20Q3I) 
models. The thick lines correspond to IRAC SI AB magnitudes, while the 
thinner lines represent K (Vega) magnitudes. The solid lines, which are 
nearly indistinguishable from one another, represent solar metallicity mod- 
els, with a Salpeter IMF. The four lines are all 13 Gyr models, with a) single 
stellar population (ssp), no dust; b) ssp, r v = 1 extinction; c) constant star 
formation rate (SFR), no dust; and d) constant SFR, t v = 1 extinction. The 
two models with dashed lines are constant SFR models, at supersolar and 
subsolar metallicity; these models provide the most different k-corrections. 
There is almost no sensitivity to population age or initial mass function (not 
shown). 

the image that are not as well centred. For each group in the sample, 
therefore, we compute how many of the spectroscopic members 
have NIR data from either set of observations. Groups are only in- 
cluded in the sample if the fraction of members with NIR data is 
at least 2/3, and for each group the members are weighted by the 
inverse of this fraction. These groups are listed in TableQ with the 
number of spectroscopically confirmed group members with either 
INGRID or Spitzer observations. There are 29 (38) groups with 
good INGRID (IRAC) coverage, and a total of 58 that have cov- 
erage with at least one of the two instruments. When a galaxy is 
observed with both instruments, we use the IRAC data since it has 
the more stable zeropoint. 



2.7 Combined infrared sample 

To present luminosity functions it will be useful to have all the 
infrared data (taken in two different filters and over a range of 
redshifts) corrected to the same rest-frame wavelength range. We 
will take the approach, described in detail below, of converting 
the Spitzer SI AB magnitudes to equivalent A"— band (Vega) at 
z — 0.4, for ease of comparison to the literature. This magni- 
tude will be denoted K^ A and is on the Vega system; note that 

KaB ~ AVcga + 2.0. 

First, k-corrections are computed to convert the INGRID and 
IRAC magnitudes separately to the corresponding observed wave- 
length range at a fiducial redshift. To keep the correction as small 
as possible, where convenient we will show the data k-corrected to 
they typical group redshift of z — 0.4, rather than z — 0. To make 
these corrections we compute the difference in K magnitude rela- 




2.5 

( S1 0.4_ K 0.4) 

Figure 4. The distribution of (SI 0,4 — K 0A ) colours for the 39 galaxies 
observed with both INGRID and IRAC, and with 0.3 < z < 0.6, K < 
17 and 51 < 20. This colour is insensitive to the nature of the stellar 
population, so we use the mean colour (2.45) as a conversion between the 
two magnitudes. The solid line is a Gaussian with this mean, and standard 
deviation 0.17 mag, a best fit to the 3-cr clipped data. 



tive to z = 0.4 for a variety of Bru zual & Charlotl 12003) galaxy 
models. The results are shown in Figure [5] and the model parame- 
ters are described in the caption. Since the NIR fluxes are relatively 
insensitive to these parameters, the choice is not critical, and to a 
good approximation we find we can use 



K(z) - K 



Sl(z) - Sf 







0.3 < z < 0.6 



-1.732 + 0.52 2<0.3 



-2.502 + 1.0 0.4<z<0.6 
-1.82 + 0.72 2 <0.4 



(5) 



(6) 



where K and 51 ' refer to the observed band at z — 0.4. 

After applying these corrections we can compare the K 0A and 
SI ' 4 magnitudes for galaxies observed with both instruments. In 
particular there are 39 such galaxies, with 0.3 < z < 0.6, K < 17 
and SI < 20 (one magnitude brighter than the photometric com- 
pleteness limit, to reduce the impact of magnitude uncertainties). 
The distribution of (SI ' -K° ) is shown for these galaxies in Fig- 
ure|4| After clipping the only 3-cr outlier, the mean and median are 
both 2.45, and the standard deviation is 0.17 mag . There is no sig- 
nificant trend of this difference with rest-frame optical colour or 
redshift. Thus we take the dispersion to represent the uncertainty 
in photometric calibration, including k-corrections and aperture ef- 
fects. We convert our SI ' 4 magnitudes to an equivalent K 0A by 



A 



sr 



2.45. 



(7) 



2 We note that this colour is also in excellent agreement with predictions of 
the same models used to calculate the k-corrections: all models considered 



predict (SI ' 4 - K 0A ) = 2.4 ± 0.1. 
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Table 3. Average group properties, divided into redshift (column 1) and velocity dispersion (a, column 2) bins. The other columns are: (3) The number of 
groups contributing to each bin; (4) the total observed, statistically weighted, rest-frame K-band luminosity within R200 and brighter than M*\ (5) The total 
rest-frame K-band luminosity including a correction for galaxies below M* , assuming a = —0.8; (6) the same as (5), but for a = —1.09; (7) the average 
dynamical mass in this bin; (8), the mass-to-light ratio, defined as column (7) divided by column (5). 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Redshift 


cr 


-/V(j TOWpS 


Lk, M k < M* 


L K ,tot (a = -0.8) 


L Kt tot (a = -1.09) 


-M200 


M200 / L>K,tot 




(kms^ 1 ) 




(10 n L K , Q ) 


(10 n £j<-,s) 


(lO 11 /^,©) 


(1O 13 M ) 


(M q /L k 'q) 


0.1-0.25 


0-250 


7 


2.03 ±0.4 


4.48 ±0.8 


6.17 ± 1.1 


1.06 


23.77 ±4.2 




250-425 


4 


3.08 ±0.7 


6.8 ± 1.6 


9.36 ± 2.2 


3.25 


47.79 ± 11.4 


0.25-0.37 


0-250 


5 


4.23 ± 1.5 


9.33 ±3.4 


12.85 ±4.7 


0.63 


6.79 ±2.5 




250-425 


7 


4.96 ± 1.5 


10.94 ±3.4 


15.07 ±4.6 


4.04 


36.91 ± 11.4 




425-700 


7 


3.41 ± 1.2 


7.52 ±2.6 


10.36 ±3.6 


18.04 


239.88 ± 84.4 


0.37-0.6 


0-250 


9 


2.97 ±1 


6.55 ±2.3 


9.02 ± 3.1 


0.39 


5.99 ±2.1 




250-425 


11 


6.87 ±0.8 


15.17 ±1.7 


20.9 ± 2.3 


3.64 


24.01 ± 2.6 




425-800 


9 


9.19 ±2.2 


20.28 ±4.9 


27.93 ±6.7 


21.69 


106.95 ± 25.8 


0.1-0.6 


0-250 


21 


2.95 ±0.6 


6.52 ± 1.4 


8.98 ± 1.9 


0.75 


11.55 ± 2.5 




250-425 


22 


5.57 ± 1 


12.3 ± 2.2 


16.95 ±3.1 


3.67 


29.81 ± 5.4 




425-800 


15 


7.1 ± 1.4 


15.68 ±3.2 


21.59 ±4.4 


20.37 


129.92 ± 26.4 



Finally, from Equation[3]we see that the rest-frame K-band magni- 
tudes are well approximated as 

tfrest = K 0A + 0.52 (8) 

= K° cq 4 + 0.52. 

This k-correction between 2 = 0.4 and z — is good to within 
about 0.05 mag, for all solar metallicity models considered. 



3 RESULTS 

3.1 Luminosity Function 

Most of our groups have less than about ten members with red- 
shifts. To construct a luminosity function it is therefore necessary 
to stack the groups 3 . When computing the luminosity functions we 
consider only galaxies within R200, which itself is estimated from 
the velocity dispersion of the group (Equation |2j- Absolute mag- 
nitudes are computed from the i\" rC st magnitudes described above, 
according to our adopted cosmological model. 

Since our data combine two spectroscopic samples, with dif- 
ferent completeness limits, we first compute the luminosity func- 
tion separately for the R21.5 and R22 samples, and then combine 
them, appropriately weighted according to the relative number of 
contributing galaxies. We have verified that the luminosity func- 
tion for either sample alone is statistically consistent with that of 
the combined data. Since the R22 sample is more complete, this 
suggests that any systematic effects resulting from undersampling 
are masked by the larger, statistical uncertainties. 

The luminosity function for all groups in the survey is shown 
in Figure|5] The number of galaxies per magnitude bin are weighted 
by the appropriate statistical weights, and divided by the number of 
groups contributing to that bin. The dotted line shows the lumi- 
nosity corresponding to the magnitude limit of the highest redshift 

3 As discussed in Section 12.31 the spectroscopic completeness is not a 
strong function of group size, which ensures that this stacked group is not 
dominated by sparsely-populated groups with large statistical weights; the 
weights are primarily a function only of apparent magnitude. 
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Figure 5. The luminosity function, per group, for all groups in the sample. 
The error bars are computed assuming Poisson statistics on the number of 
(unweighted) galaxies in each bin. The vertical, dotted line represents the 
brightest magnitude limit of all the groups in the sample; to the right of 
that line the data are still statistically complete, but there are progressively 
fewer contributing groups. The x-axis is the absolute rest-frame magnitude 
in the K filter. The curved, solid (dashed) line is a Schechter function with 
M* = -24.38 and a = -1.09(-0.8). 

group in the sample; at fainter luminosities, fewer groups contribute 
to each bin. 

We will not attempt to derive Schechter function parameters 
from our data, which are not much deeper than M* and include 
few spectroscopic group members when the groups are binned by 
velocity dispersion (see below). However, it is necessary to as- 
sume some shape to correct the total group luminosity for galax- 
ies below the magnitude limit. iDrorv et alj 120031) have measured 
the field K-band luminosity functions out to z — 1, and found 
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Figure 6. The cumulative luminosity per group, binned by group redshift and velocity dispersion as shown. Note that in this representation, the observational 
eiTor bars (not shown) are not independent, and are largest at the brightest magnitudes. The number of groups contributing to each redshift and velocity 
dispersion bin is also shown in each panel. The x-axis is the rest-frame K-band absolute magnitude. The curved, solid (dotted) line is a Schechter function with 
MZ = —24.38 and a = — 1.09(— 0.8). Both functions are normalized so that the total luminosity in galaxies brighter than AfJ matches the observations. 



(for h = 0.7) Ml = -24.16 - 0.53z; at z = 0.4, therefore, 
M'l — —24.38. We will adopt this value throughout the paper, 
although we might expect M* t o vary by ~ 0.2 m ag over the 
full redshift range of our sample. Dr orv et alJ {2003) cannot con- 
strain the faint end slo pe, and s o adopt a = —1.09 measured from 
2MASS iKochanek et alfcOOll) . However, the slope may be some- 
what shallower in dense environments, more like a — —0.8 (e.g. 
iBalosh et alj200lURamella et alfeool) . so we will show both. The 
Schechter function using these parameters is shown in Figure [5] 
normalized to match the total luminosity of galaxies brighter than 
MJf ; both provide a reasonable description of the data, within the 
uncertainties. 



3.2 Total group luminosities 

In Figure [6] we show the cumulative luminosity distribution in 
groups, divided into redshift and group velocity dispersion bins. 
This alternative representation of the luminosity function shows the 
integrated K-band luminosity per group, brighter than a given mag- 



nitude. Models are shown based on the Schechter functions pre- 
sented in the previous section (Ml — —24.38 and a = —1.09 
or a = —0.8), always renormalized to match the total luminosity 
of all galaxies brighter than Ml, There is no dramatic change in 
the shape of the distribution with redshift or velocity dispersion, al- 
though an evolut ion of Ml o f ~ 0.2 mag over this redshift range, as 
found bv lDrorvetaTi r2003). would still be consistent with the data. 
For all bins, the data extend to Ail or deeper. To calculate the total 
luminosity, we will take the integral of the Schechter function fit; 
because of the way this is normalized, this is equivalent to adding 
up all the observed light brighter than Ml (including appropriate 
statistical weights), and using the Schechter function to extrapolate 
to zero luminosity. With this method of normalization the total lu- 
minosity is insensitive to the choice of A/*, and choosing a weakly 
evolving value would have no influence on our results. However, 
the integral does depend on the assumed value of a, so we calcu- 
late results for both a = —1.09 and a = —0.8. These results, as 
well as the total observed luminosity per group (uncorrected for the 
magnitude limit), are tabulated in Table|5| 
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Figure 7. The total near-infrared luminosity of galaxy groups is shown as 
a function of dynamical mass, as estimated from the group velocity dis- 
persion. The small, open symbols are low-redshift data from lRamella et at] 
i 2004). The filled symbols are our groups, in three bins of velocity disper- 
sion. The dynamical mass plotted is the weighted average mass of all groups 
in the bin. The error bars on this quantity show the minimum and maximum 
mass of the groups in this bin. The total stellar mass includes an extrap- 
olation to zero luminosity, assuming a = —0.8. The triangles include a 
small correction for passive evolution, based on a single stellar population 
model from lBruzual & CharloJ |2003|) . and for hierarchical growth (in both 
the stellar and dark mass) from the Millennium simulation. The solid lines 
indicate lines of constant A/200 /Lk ratio. 

The uncertainty in the total luminosity assumes the un- 
weighted number of galaxies per magnitude bin is described by 
a Poisson distribution. This accounts for the fact that most of the 
luminosity comes from the few brightest galaxies, which are gener- 
ally the least numerous and therefore subject to the largest statisti- 
cal fluctuation within a population of similar groups. It also accom- 
modates the possibility that occasionally the brightest group galaxy 
may not have a spectroscopic redshift (as long as this galaxy isn't 
systematically overlooked in every group, and is not much brighter 
than the brightest galaxy for which we do have a redshift). Fur- 
thermore, the uncertainty scales with the statistical weight on the 
galaxy, so any incompleteness at the bright end of the luminosity 
function will be reflected in correspondingly larger error bars. In 
fact, our spectroscopic completeness is very high for the brightest 
galaxies, especially for the R22 sample, and therefore we are un- 
likely to have missed the brightest member galaxy in most groups. 

3.3 Mass-to-light ratios 

We are now in a position to compare the total NIR luminosity with 
the dynamical mass of the galaxy groups. The latter quantity is 
compu ted from the v elocity dispersion (Equation in the same 
way as lRamella et alj 120041) . who performed a very similar exper- 
iment at 2 ~ (see below). The large, filled circles in Figure Q 
represent our data, shown in three bins of velocity dispersion, but 
summed over all redshifts, as in the last three rows of Table [3] 
The filled triangles demonstrate the effect of applying two, rela- 
tively small, evolutionary corrections to our data. The first is an es- 



timated correction for passive luminosity evolution, based o n a sin- 
gle stellar population, dust-free iBruzual & Charlol 120031) model, 
that amounts to a ~ 30 per cent reduction in luminosity, at fixed 
mass. The second correction aims to account for the hierarchical 
mass growt h to z = 0. Us ing the results of the Millennium simula- 
tion ISpringel et all2005l) . we find that haloes with M > 1O 13 M 
grow by ~ 43 per cent between z = 0.4 and z — 0, independent 
of initial mass. We assume that this growth occurs in both the lu- 
minosity and mass of the system, so the points move up and to the 
right, along lines of constant mass-to-light ratio. The net result of 
these two corrections is that the A/200/ 1 Lk of our 2 = 0.4 groups 
are best compared to the M200 /Lk of 2 = systems that are ~ 
40 per cent more massive, as these represent, statistically, the most 
likely descendents. 

T he small, open s ymbols in Figure^represent the low-redshift 
data of iRamella etafl 120041) . This is a sample of groups selected 
from a complete spectroscopic survey at 2 « 0, followed up 
with subsequent, dee per spec troscopy and cross-correlated with 
the 2MASS I Jarre tt et al jiooTJ) to obtain near-infrared magnitudes. 
Thus the group selection is similar to ours, in that it is based on 
redshift-space clustering rather than X-ray or radio emission. There 
are naturally differences in our procedure for identifying groups 
and measuring their velocity dispersions, in part due to the fact 
that our galaxies are drawn from an incomplete redshift survey and, 
therefore, must be statistically weighted. However, due to the small 
number of members per group, our velocity dispersions are domi- 
nated by statistical uncertainties, rather than any systematic effects 
related to exactly which gal axies w ere included in the computation. 
For example. lWilman et alJi2005al) find that the deeper, more com- 
plete spectroscopic sampling from Magellan certainly improves the 
statistical uncertainties of our velocity dispersions, but does not re- 
veal a systematic bias relative to measurements made from the orig- 
inal survey. Since we use the same equations to relate velocity dis- 
persion to dynamical mass as lRamellaetai]l2004l) . it is, therefore, 
possible to compare their data with ours to look for evolutionary 
trends. 

There are two remarkable points about the result in Fig- 
ure |7| The first is tha t our data lie within the scatter of the low- 
redshift data o f lRamella et all f 2004). There has, therefore, been no 
strong evolution in the stellar mass fraction of galaxy groups since 
2 = 0.4. Indeed, as is evident from the data in Table|3| we see no 
evidence for significant evolution within our sample itself, which 
spans the redshift range 0.1 < 2 < 0.6. Given the statistical and 
systematic uncertainties we can rule out a pure luminosity evolu- 
tion of a factor ~ 2 or greater since 2 ~ 0.4. Also interesting is the 
suggestion that the M200/ Lk ratio at 2 = 0.4 depends on system 
mass. The lowest mass groups in our sample have M200 / Lk^ 10, 
while the most massive systems are a factor ~ 10 larger. Again the 
large uncertainties preclude strong conclusions, but this is similar 
to the trend observed at low-redshift iRamella et alj2004LlLin et"aT] 
2004). We will discuss this further in Section|4| 



4 DISCUSSION 
4.1 Lensing masses 

Dynamical mass estimates based on velocity dispersion can be 
problematic for groups, both because of the small number of galax- 
ies involved, and because of the difficulty of testing the neces- 
sary assumptions required to link velocity dispersion to mass. It 
is unclear, for example, how close these groups are to dynamical 
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Table 4. Similar to Table l3l but where we have included an estimate of the average lensing mass, <r lcns , in column (3), for each range of re dshift-derived 
velocity dispersions (column 1). The lensing masses are derived from a fit between the lensing and dynamical velocity dispersion estimates in IParked |2005|> . 
and are used to calculate the cluster masses, i?200 . ar >d luminosity within -R200 ■ Their uncertainties are the statistical uncertainties given in Parker 1 200!$). The 
other columns are as in TablelSI The uncertainty in M200 / Lk (column 8) ignores any uncertainty in the mass, which would be correlated with the uncertainty 
in L K . 
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equilibrium. Weak lensing provides a promising alternative, as this 
method measures the mass independently of the dynamical state of 
the group. Of course, this method has its own drawbacks, which 
are mainly the weakness of the signal, the sensitivity to choice of 
group centre, and the sensitivity to all mass projected along the line 
of sight between the observer and the lensed galaxy. Nonetheless, 
as the sources of systematic uncertainty are quite different from 
those associated with the velocity dispersions, they provide a use- 
ful check on our results. 

Weak lensin g masses have been m eas ured for the C NOC2 
group sample by Hoekstra et al. (2001) and Parker et al. 12005). 
However, due to the weakness of the lensing signal, this can only 
be done for ensembles, and Parked i2005l) provide masses for three 
classes of groups, split by the velocity dispersion (as measured by 
ICarlberg et a7]|2001l) at a = 190 km s" 1 and a = 500 km s" 1 . 
We make a linear fit to relate the lensing dispersions <7i cns (the 
velocity dispersion associated with an isothermal sphere of the 
measured mass) to the average dynamical values a^ yn , an d find 4 
ficns = 113km s _1 + 0.49<7d yll . Using this fit as an estimate of the 
lensing mass for each individual group, we recompute -R200 and 
M200 and repeat the entire analysis to measure the corresponding 
total NIR luminoisty. This includes recomputing the NIR complete- 
ness of each group, since this depends on R200, as well as the total 
luminosity within ifeoo- The results are shown in Table|4] 

The main difference from the results in g|3T2lis that the masses 
of the groups with the highest velocity dispersion are lower than es- 
timated using Equation^ This leads to a lower value of A/200 /Lk 
(despite the corresponding reduction in R200) and reduces the in- 
crease of M200 / Lk with system mass to a factor ~ 3 between the 
lowest- and highest-mass groups in our sample. 
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Figure 8. The cumulative stellar mass function for the whole group sample 
(top left panel) and for samples subdivided into bins of velocity dispersion 
(remaining panels, as labeled). The stellar masses are computed from the 
K res t values, assuming the simple, colour-dependent M ste n ar /Ljf model 
described in the text. Each panel shows the number of groups contributing 
to each bin, and the value of M* assumed for the Schechter functions, plot- 
ted as curved lines. The solid line assumes a faint end slope of a = —1.09, 
while the dotted line assumes a = —0.8. The Schechter functions are nor- 
malized to match the observed total stellar mass in galaxies with M > M* . 



4.2 Stellar mass-to-light ratios 

While the near-infrared luminosities are closely related to stel- 
lar mass, there is still substantial variation in the stellar mass- 
to-lig ht ratio (M s t e ii a r I 'L k) with galaxy colour. We use the 
iBruzual & CharloU 120031) models to compute M s t e iiar / Lk as a 
f unction of galaxy colour, for a range of parameters. Assuming 
a IChabrieil 120031) initial mass function, we model blue galaxies 
(B — V < 0.4) as young galaxies with constant star formation 
and dust extinction of t v = 1 mag; this gives i\/ stc n al ./ Lk— 0.2. 
For red galaxies (B — V > 1), a 11.7 Gyr old, dust-free single 

4 To be fully c onsist ent, we use the original measurements o f oyivn from 
Carlbe rg et alj 1200 ll) to bin the groups in the same way as IParker et alj 
1 2005), but then use our updated values to compute the average o"<j yn within 
each bin. 



stellar population model gives M a t B iiar/£if = 0.7. For intermedi- 
ate colours, we simply adopt a linear interpolation between these 
values. 

The cumulative stellar mass functions are shown in Figure [8] 
for the whole group sample and for subsamples divided by ve- 
locity dispersion. The mass functions can be reasonably approxi- 
mated by the Schechter functions shown, which have characteristic 
mass scales ranging from 10 10 ' 9 Mq for the lowest-mass groups, 
and lO 11,1 ^/© for the highest-mass group s. This trend for M * to 
increase with density is well-known (e.g. Baldrv et al. 200^). For 
the entire group sample, the ratio of dynamical mass to total stellar 
mass is 168±25. In Figure|9|we show the dynamical-to-stellar mass 
ratio as a function of velocity dispersion, for both the a— based and 
weak lensing-based mass (and R200) estimates. We also show the 



The stellar mass content of distant galaxy groups 1 1 




> 
i 

m 



X. X • 



x x x x 




• ■ x X 



« XX 



Groups 



• x x ■ x x 

: V xx |x v** ** x 

.•UxX* xxx xx x 



a(km/ s) 

Figure 9. The dynamical-to-stellar mass ratio as a function of velocity 
dispersion. The solid points are computed assuming dynamical masses 
and R200 derived from the velocity distribution, while the open sym- 
bols use the weak-lensing derived masses. The circles are based on our 
simple colour-dependent M ste n ar /Lx model, while the triangles assume 
M s toiiar/iif = 0.7 for all galaxies. The three velocity dispersion bins are 
the same as in Figure [5] and the points are shown in the middle of each 
bin with the range indicated by the error bars on one set of points only, for 
clarity. 
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Figure 10. The correlation between (B-V) rest-frame colour and stellar 
mass, for the group (bottom panel) and field (top panel) samples. Both sam- 
ples are restricted to the redshift range 0.12 < 2 < 0.55, where there is no 
colour incompleteness. Only galaxies brighter than K = 17.2 (correspond- 
ing to R = 21.5) are shown (crosses), and group galaxies are restricted to 
those within -R200- Filled circles represent galaxies with strong emission 
lines, Wo([OlI])> IOA. The few small points represent galaxies for which 
no Wo ([Oil]) measurement is available. 



effect 5 of using a constant Miliar / Lk = 0. 7, as expected for early 
type galaxi es, and typ ical of the Universal average at low redshift 
(e.g. lCole et all200lf) . The spread in the points at a given velocity 
dispersion therefore gives some indication of the systematic un- 
certainty on the measured M200 /-^stellar- As already noted, there 
is some discrepancy between the lensing-based and redshift-based 
results for the systems with the highest velocity dispersions, but 
overall the data show that M200 /-Msteliar increases from g 50 in 
the poorest systems to Z 100 in the massive groups. This could be 
interpreted as a factor ~ 2 difference in star formation efficiency 
between these haloes. 



4.3 Colours and emission line strengths 

In IWilman etafl <2005ah . we found that the fraction of galaxies 
with [Oil] emission increased with decreasing, rest-frame Bj-band 
luminosity, for both the group and field populations in the CNOC2 
(R22) galaxy sample. We also found that, at fixed luminosity, emis- 
sion line gal axies were more comm on in the field population. Using 
toy models (Wilman et al. 2005b) we claimed that this difference 
was due to a recent truncation of star formation in group galaxies, 
qualitatively similar to wha t is observed in more massive c lusters 
at this redshift (e.g. |BaloghetalJll99l. l200l iNakata et alJl2005t 
IPoggianti et alj Eo06). However, that analysis was complicated by 



5 The difference between these two measurements is not simply related 
to the difference in average A^tellat/^Ki because the fixed value of M* 
means the normalization of the Schechter function, and hence the extrapo- 
lation to Lk = 0, also differs in the two cases. 



the fact that the rest-frame Bj-band luminosity is sensitive to re- 
cent star formation. Here we can re-examine the correlation as a 
function of stellar mass. 

Figure [Tol shows the rest-frame (B-V) colour as a function of 
stellar mass. We restrict the redshift range to 0.12 < z < 0.55; 
outside this range, colour-dependent incompleteness becomes im- 
portant because few identifiable absorption features of red galax- 
ies lie within the band-limiting filter used for the original CNOC2 
spectroscopy <Yeeetall2 000). Within this range, all galaxies with 
redshifts are significantly detected in all filters, so there is no colour 
bias. Galaxies with strong emission lines (Wo([Oll])> 10A) arc 
shown as filled symbols. 

Both the group and field populations exhibit a decrease in the 
fraction of blue, emission line galaxies as stellar mass increases. 
However, at fixed stellar mass, such late-type galaxies are less com- 
mon in the group environment, compared with the field. This is il- 
lustrated more clearly in Figure fTTI where we show the fraction of 
galaxies with strong emission lines (Wo([Oll])> 10A) as a func- 
tion of stellar mass, for the field and group samples. Thus we con- 
firm our conclusions in Wi lman et alj l2005bl) . that star formation 
in group environments was terminated earlier than in the lower- 
density field at z ~ 0.4. 

4.4 Interpretation 

The dominance of old stellar populations, together with the weakly- 
evolving mass-to-light ratio, implies that most of the stellar mass in 
groups must have already been in place by z ~ 0.4. Thus, we see 
no evidence of accelerated star formation in groups, despite the fact 
that the low velocity dispersions are expected to lead to enhanced 
dynamical interactions between group members. This is especially 
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Figure 11. The fraction of galaxies with strong [Oil] emission lines is 
shown as a function of stellar mass, for the field and group sample at 
0.12 < z < 0.55. Error bars are 1-cr jackknife resampling estimates. 



surprising since our sample is at an epoch when the global star 
formatio n rate was several times higher than it is at present (e.g. 
lHopkinslE 004). Unlike with group samples selected from their X- 
ray emission (e.g. iMulchaev et aljl200rJ) . we have little reason to 
expect our redshift-selected sample to be biased toward especially 
relaxed or evolved groups, and therefore conclude that older stellar 
populations are characteristic of groups at z ~ 0.4 in general. 

Recent numerical simulations suggest that in overdense envi- 
ronments, dark matter haloes assembled their mass more rapidly, 
and at higher redshift, than haloes of the same mass in low-density 
environments iGao et all2005h . and this might be sufficient to ex- 
plain the more evolved populations of galaxy groups and clusters 
without appe aling to local interaction s between galaxies and their 
environment l Maulbetsch et al. 2006). Determining whether or not 
this is a viable explanation for the observations presented here 
awaits a quantitative comparison with theoretical models. 



5 CONCLUSIONS 

We have presented near-infrared observations of 58 redshift- 
selected galaxy groups at 0.1 < z < 0.6, obtained from the WHT 
and Spitzer space telescope. This affords us the first opportunity to 
obtain reliable stellar masses for galaxies in groups at this redshift. 
From these data we draw the following conclusions: 

• There is evidence that the highest-mass groups in our sample 
(<r > 500 km s _1 ) have M200 / Lk that are a factor Z 3 larger than 
for the lowest-mass groups (a < 200 km s _1 ). This trend is present 
whether we use dynamically-estimated masses (from the velocity 
dispersions) or masses measured from a weak lensing analysis. 

• Our best estimate of the group stellar mass fraction decreases 
from ~ 2 per cent at a < 200 km s" 1 , to < 1 per cent in the most 
massive groups with a > 500 km s -1 . 

• When comparing groups at z ~ 0.4 to their statistically most- 
likely descendants at z = 0, we find no evidence for strong evo- 



lution in M200 / Lk beyond that expected for a passively evolving 
population. 

• Group galaxies have older stellar populations, as measured 
by their ^([Oll]) emission or rest-frame (B-V) colour, than field 
galaxies of the same stellar mass. 

We have demonstrated that galaxies in groups are predomi- 
nantly old, passively evolving systems, even at z ~ 0.4. This is 
something of a surprise, given that effects such as tidal interactions 
and mergers (which are expected to induce star formation) should 
be more common, while more dramatic effects like ram-pressure 
stripping, often thought to play a role in terminating star formation 
within more massive clusters, should be negligible. Our results im- 
ply that these effects have not had a large influence on the properties 
of galaxies in groups since z ~ 0.4. 
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